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Abstract 

The linear complete differential resultant of a finite set of linear or- 
dinary differential polynomials is defined. We study the computation 
by linear complete differential resultants of the implicit equation of a 
system of n linear differential polynomial parametric equations in n — 1 
differential parameters. We give necessary conditions to ensure proper- 
ness of the system of differential polynomial parametric equations. 



1 Introduction 

The implicitization problem of unirational algebraic varieties has been widely 
studied and the results on the computation of the implicit equation of a sys- 
tem of algebraic rational parametric equations by algebraic resultants are well 
known [9], [10]. This work was motivated by the paper of X.S. Gao [11], where 
computational issues, related with the implicitization problem of differential 
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rational parametric equations, are treated by characteristic set methods, mak- 
ing use of the differential algebra theory developed by Ritt [ISj and Kolchin 
[12]. The paper by X.S. Gao [11], establishes the basis ground for the general- 
ization to the differential case of the results in algebraic geometry on implicit 
and parametric representations of unirational varieties, conversion algorithms, 
etc (see, for instance, [5], [TU], [IB])- 

We explore the first steps of the generalization to the differential case of 
the results in algebraic geometry on implicit representations of unirational va- 
rieties. To be more precise, we are interested in finding a differential resultant 
that would solve the differential implicitization problem. The implicitization 
problem of differential rational parametric equations is a differential elimina- 
tion problem. There is a wide range of applications of the differential elimina- 
tion method to computer algebra and applied mathematics. For a survey on 
differential elimination techniques and their application to biological modeling 
we refer to [3] . 

We defined the implicit equation of a system of n differential rational para- 
metric equations in n — 1 differential parameters in [T7]. In this paper, we 
study the computation by differential resultants of the implicit equation of a 
system P(X, U) of n linear differential polynomial parametric equations (lin- 
ear DPPEs) in n — 1 differential parameters mi, . . . ,Mn-i (we give a precise 
statement of the problem in Section [2]), 



The differential resultant of a set of ordinary differential polynomials was 
introduced and studied by G. Carra'Ferro in [5] for two differential polynomials 
and in [6] for a finite number n of differential polynomials in n — 1 differen- 
tial variables. The generalized differential resultant was also defined by G. 
Carra'Ferro in [7] for n + s differential polynomials in n — 1 differential vari- 
ables, s > 0. Previous definitions of the differential resultant for differential 
operators are due to L.M. Berkovich-V.G. Tsirulik [2] and M. Chardin [8]. 

Let us consider the linear ordinary differential polynomials Fi{X,U) = 
Xi — Pi{U) of order Oj, i = 1, . . . , n. As we will explain in Section 13.11 the 
differential resultant c}Res(Fi, . . . , F„) is the determinant of the L x L matrix 
M{L), where L = ^"=i((^fc=i o^) — Oi + 1). In Section HI we prove that 
if (9Res(Fi, . . . , F„) 7^ then the imphcit equation of the system V{X, U) is 
given by the differential polynomial c)Res(-Fi, . . . ,F„)(X). Then we analyze 
some of the reasons for 9Res(Fi, . . . , F„) = 0. 

The differential resultant 9Res(Fi, . . . , F„) is the Macaulay's algebraic re- 
sultant of a differential polynomial set PS with L elements. One reason 
for 9Res(Fi, . . . , Fn) = is the following: if the polynomials in PS are not 



V{X,U) 




PniU) 
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complete in all its variables then the resultant will be zero. If for some 
j G {1, ... ,n — 1} the differential variable Uj has order less than Oj in Fi 
for all i G {1, . . . ,n} then the matrix M{L) has one or more columns of zeros. 
Let 7 be the number of zero columns in M{L) due to this reason. Let us see 
an illustrating example. 

Example 1.1. Let us consider the following system of linear DPPEs 

Xi = Ui + U2 + U2l 
X2 = tUu + U22 
X3 = Ui + U21 

where Ujk = d^Uj/dt^, j = 1,2 and k ^ N. The differential resultant of 
Fi{X, U) = X1-U1-U2-U21, F2{X, U) = X2-tuii-U22, Fsi^X, U) = X3-U1-U21 
is zero because the order ofui in every polynomial is less than its order, 7 = 1. 
The implicit equation of this systems is {t — l)xi2 — tx^i — {t — l)x32 + X2, 
Xif^ = d^Xi/dt^, z = 1, 2, 3, A; G N. 

We define in Section [5] the linear complete differential resultant of a finite 
set of linear ordinary differential polynomials, generalizing Carra'Ferro's dif- 
ferential resultant in the linear case. We prove the next theorem in Section 

El 

Theorem. Given a system V{X, U) of differential polynomial parametric 
equations. If the linear complete differential resultant dCR.es{Fi, . . . , F^) 7^ 
then 9CRes(-Fi, . . . , F„)(X) = is the implicit equation ofV{X, U). 

In Section [TJ we study whether the system U) of linear DPPEs is 

proper, problem closely related to the existence of inversion maps. If the sys- 
tem U) is not proper we prove that 9CRes(Fi, . . . , F„) = 0. The homo- 
geneous part of the linear differential polynomials Fi, . . . , can be written 
in terms of differential operators. We obtain necessary conditions on these 
differential operators so that c}CRes(Fi, . . . , F^) 7^ 0. 

Computations throughout this paper were carried out with our Maple im- 
plementation of functions to compute the differential resultant by CarraTerro 
in [6] and the linear complete differential resultant defined in Section El of this 
paper. Our implementation is available at [TB^. 

The paper is organized as follows. In Section [51 we introduce the main no- 
tions and notation. Next we review the definition of the differential resultant 
defined by Carra'Ferro in Section [31 and we define the differential homogeneous 
resultant. In Section [Hwe explain the computation of the implicit equation 
by CarraTerro's differential resultant. We define the linear complete differ- 
ential resultant in Section [SI In Section [HI we give our main results on the 
implicitization of linear DPPEs by linear complete differential resultants. Our 
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results on properness appear in Section [71 We finish the paper discussing the 
special cases n = 2 and n = 3 in Section [HI 



2 Basic Notions and Notation 

In this section, we introduce the basic notions related to the problem we deal 
with, as well as the notation and terminology used throughout the paper. For 
further concepts and results on differential algebra we refer to [12] and [T5] . 

Let IK be an ordinary differential field with derivation d, ( e.g. Q(t), d = 
^). Let X = {xi,...,Xn} and U = {ui, . . . , be sets of differential 

indeterminates over K. Let E be a universal extension field of K containing the 
set of indeterminates U. Let No = {0, 1, 2, . . . , n, . . .}. For e Mo we denote 
by Xik the k-th derivative of Xj. Given a set Y of differential indeterminates 
over IK we denote by {Y} the set of derivatives of the elements of Y, {Y} = 
{d^y \ y eY, k E Nq}, and by IK{X} the ring of differential polynomials in 
the differential indeterminates xi, . . . , x^, that is 



Analogously for K{[/}. 

A system of differential rational parametric equations (system of DRPEs) is 
a system of the form 



where Pi, . . . , P„, Qi, . . . , Q„ G IK{[/}, Qi are non-zero, and not all Pj, Qi G IK, 
i = 1, ... ,72. We call the indeterminates U a set of parameters of Q{X, U); 
they are not necessarily independent. When all differential polynomials Pj, Qi 
are of degree at most 1, we say that ([1]) is a linear system. Moreover, if 
all Qi e K, we say that ([1]) is a system of differential polynomial parametric 
equations (system of DPPEs). Associated with the system ([1]) we consider the 
differential ideal (see [TT], Section 3) 



and we call it the implicit ideal of the system ([T]). By [TT], Lemma 3.1, the 
implicit ideal ID is a differential prime ideal. Moreover, given a characteristic 
set C of ID then n— | C | is the (differential) dimension of ID, by abuse of 
notation, we will also speak about the dimension of a DRPEs system mean- 
ing the dimension of its implicit ideal. The parameters U are independent if 



K{X} = K[xik M = 1, . . . , 



n, k E No]. 




(1) 



ID = {/ G K{X} I f{P,{U)/Q^{U), Pn{U)/Q^{U)) = 0} 



dim(Q(X,?7)) = \U\ (see [H], Section 5). 
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The differential variety defined by 

Zero(ID) = {r/ G E" | V/ G ID, f{r]) = 0} 

is called the implicit variety of Q{X,U). If the parameters of Q{X,U) are 
not independent and K = Q(t), there exists a set of new DRPEs with the 
same implicit variety as Q{X,U) but with independent parameters, (see [11], 
Theorem 5.1). 

If dim(ID) = n — 1, then C = {A[X)} for some irreducible differential 
polynomial A G ]K{X}. The polynomial A is called a characteristic polynomial 
of ID. Furthermore, if B is another characteristic polynomial of ID then A = 
bB with 6 G K. 

We introduced the notion of implicit equation in [17] and we include it here 
for completion. 

Definition 2.1. The implicit equation of a{n — 1)- dimensional system of DR- 
PEs, in n differential indeterminates X = {xi, . . . ,Xn}, is defined as the equa- 
tion A{X) = 0, where A is any characteristic polynomial of the implicit ideal 
ID of the system. 

In the sequel, we consider the linear system of DPPEs 

f xi = Pi(t/) 
P(X, U) = l : (2) 

[ x„ = P^{U) 

So, Pi, . . . , P„ G ]K{f/} with degree at most 1, and not all Pj G K, z = 1, . . . , n. 

Let K.[d] be the ring of differential operators with coefficients in K. For 
i = 1, . . . n and j = 1, ... n — 1, there exist differential operators Cij G K[c?] 
and constants G K such that 

n-1 

Pi{U) = at -^^Cijiuj). 

i=i 

We define the differential polynomials 

n-1 

T,(X) =Xi- ai, Hi{U) = Y^^iA^j), and Fi{X, U) = T,(X) + Hi{U). 

Given P G K{X U t/} and ?/ G X U t/, we denote by ord(P, y) the order of P 
in the variable y. If P does not have a term in y then we define ord(P, y) = —1. 
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Remark 2.2. To ensure that the number of parameters is n—1, we assume that 
for each j G 1} there exists i G {1, . . . , n} such that the differential 

operator Cij ^ 0. That is, for each j G {1, . . . ,n — l} there exists i G {1, . . . , n} 
such that ord(-Fj, uj) > 0. 

In this situation, the problem we deal with in this paper is: given a system 
V{X, U) of linear DPPEs, to compute an implicit equation using differential 
resultants. In [11], algorithmic methods for solving this problem in a more 
general case are presented in the language of characteristic sets. Our candidate 
to be the imphcit equation of V{X, U) is the differential resultant of Fi, . . . F„ 
that we describe in the next section. 

3 Differential Resultants 

Let © be a differential integral domain, and let fi G D{f/} be an ordi- 
nary differential polynomial of order Oj, z = 1, . . . , ra. A differential resultant 
9Res(/i, . . . , fn) of n ordinary differential polynomials /i, ...,/„ in n — 1 dif- 
ferential variables Ui, . . . , Un-i was introduced by Carra'Ferro in [B]. Such a 
notion coincides with the Macaulay's algebraic resultant [13] of the differential 
polynomial set 

n 

PS(/i, ...,/„):= {S^-"'/., . . . , 9/,, /, I z = 1, . . . , n, where iV = o,}, 

i=l 

Now, let hi G ro{?7} be a ordinary differential homogeneous polynomial 
of order Oj, i = l,...,n. We define the differential homogenous resultant 
dRes^{hi, . . . , hn) of n ordinary homogeneous differential polynomials hi, . . . ,hr. 
in n — 1 variables as the Macaulay's algebraic resultant of the differential poly- 
nomial set 

n 

PS^ (/ii, ...,hn):= {d^-°^-^hu . . . , dhu /li I z = 1, . . . , n, where = ^ a}. 

i=l 

A differential homogeneous resultant was defined also by Carra'Ferro in [5] 
for n = 2. In addition, when the homogeneous polynomials have degree one 
and n = 2 the differential homogeneous resultant coincides with the differential 
resultant of two differential operators defined by Berkovitch-Tsirulik in [2] and 
studied also by Chardin in [S]. 

We implemented in Maple the differential resultant matrices defined by 
Carra'Ferro in [B], our implementation is available at [IB]. This is the tool 
used to perform our experiments and in particular the computations in the 
examples of this paper. 
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Differential resultants are Macaulay's algebraic resultants therefore with 
some previous computations, the implementation of the Macaulay's algebraic 
resultant ( available at [H]) could be also used to compute differential resul- 
tants. 

Example 3.1. Let K = C{t) and d = The differential resultant c?Res(/, g) 
of the differential polynomials f{ui) = t — Auu"^ — iuf — tunUi — 5uii — iui and 
g{ui) = t — ui — Suii in K{ui} ( uu = d^Ui/df) is the Macaulay's algebraic 
resultant of the set PS{f,g) = {df,f,dg,g}. By ^B], Definition 10, we can 
compute dlles{f,g) as the quotient of two determinants. The numerator is the 
determinant of a matrix of order 20, the number of monomials below, that we 
call M(20) and the denominator is the determinant of a submatrix of M(20) 
that we call A. 

The order of the differential polynomials in PS{f,g) is less than or equal 
to two. The rows of the matrix M(20) are the coefficients of the following 
differential polynomials: 

rows 1... A {uudf,uidf,ui2df,df} 

rows 5... 8 {nil/, ui/, U12/, /} 

rows 9 ... 16 {uuUidg, uuUi2dg, uiUi2dg, uu'^dg, 

uiidg,uidg,ui2dg,dg} 
rows 17... 20 {uuuig,uug,uig,g}. 

The coefficients are written in decreasing order using first the degree and 
then the lexicographic order with U12 < Ui < Uu, that is, the columns of the 
matrix M(20) are indexed by the monomials in the list 

Mll^, Mll^Ml, Mll^Ml2, MllMi, M11M1M12, MllMl2^, u\, u\ui2, MlMl2^, ^12^, Mll^, Mll^l, 
M11M12, Ml, M1M12, Ml2^ Mil, Ml, M12, 1- 

The matrix A is the submatrix of M {20) obtained by removing rows 1, 2, 5, 
6, 10, 11, 12 and columns 1, 2, 4, 6, 7, 9, 10 o/M(20). Finally, the differential 
resultant is the quotient dKes{f,g) = det(M(20))/ det(/l) where 

det(M(20)) = 1890501 12 1 + 648075168^2 + 274613328 + 1039857264 
- 240663312 - 108661824 + 1932336 + 3114288 - 174960 
det(A) = -13608 1^ - 64476 - 48600 + 8748 1^ - 3888 1^. 



3.1 Differential resultant of Fi, . . . , F„ 

We consider now the polynomials Fi and Hi, introduced in the Section [2], 
and we set © = K{X}. In this section, we give details on the computa- 
tion of 9Res(Fi, . . . , Fn) and dRes^{Hi, . . . , i/„), since they will be important 
tools in this paper. We think of Fi, . . . , F„ as polynomials in the n — 1 vari- 
ables ui, . . . ,Un-i and coefficients in the differential domain D; recall that 
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Fi, . . . , F„ are of orders Oj and degree one. We review below the computation 
of 9Res(Fi, . . . , Fn) by means of determinants as in [B]. 

We define rankings on the sets of variables X and U so that the matrix we 
use to compute the differential resultant 9Res(Fi, . . . , F„) equals the one used 
in [6]. 

• The order x„ < ■ ■ ■ < Xi induces a ranking on X (i.e. an order on {X}) 
as follows (see [12], page 75): x < dx and x < x* ^ d^x < d^*x*, for all 
x,x* e X, k, k* e No- 

• The order ui < ■ ■ ■ < Un-i induces an orderly ranking on U as follows (see 
[12], page 75): u < du, u < u* ^ du < du* and k < k* ^ d'^u < d^* u*, 
for all M, u* e U, k, k* E Nq. We set 1 < ui. 

We call TZ the ranking on X U U that eliminates X with respect to U, 
that is d'^x > d''*u, for all x E X, u E U and k,k* E Nq. Now, the set 
PS = PS(Fi, . . . ,Fn) is ordered by TZ. Note that, because of the particular 
structure of Fi, one has that: 

Fn<dFn<---< d^-^-Fn < ■ ■ ■ < F2 < dF2 < ■ ■ ■ < Fi < ■ ■ ■ < d^-^'Fi. 

That is, PS is a chain (see [15], page 3) of differential polynomials {Gi, . . . , Gl,} 
with L = X;r=i N -Oi + l = {n- l)N + n; recall that = Y^^^^i Oi- 

Then, let M{L) be the LxL matrix whose k-th row contains the coefficients 
of the {L — k + l)-th polynomial in PS, as a polynomial in D{t/}, and where the 
coefficients are written in decreasing order with respect to the orderly ranking 
on U. Hence, M{L) is a matrix over KjX} that we call the differential resultant 
matrix of Fi, . . . , Fn- In this situation: 

<9Res(Fi, ...,Fn) = det(M(L)). 

Analogously, we can use determinants to compute dRes^{Hi, . . . , Hn); re- 
call that the homogeneous differential polynomials Hi E K{U} are of orders 
and degree one. Let = L — n and consider PS'^ = PS'^(ifi, . . . , if„) as the 
polynomial set obtained from PS by subtracting from the chosen polynomials 
its monomial in D (i.e. Xi — ai), therefore maintaining in PS'* the ordering 
estabhshed in PS. Let M(L^) be the L'^ x L'* matrix whose {L^ - k + l)-th 
row contains the coefficients of the k-t\i polynomial in PS'*, as a polynomial in 
D{f/}, and where the coefficients are written in decreasing order with respect 
to the orderly ranking on U . Hence, M{L^) is a matrix over IK that we call 
the differential homogeneous resultant matrix of Hi, . . . , Hn- In this situation: 

dRes''{Hi, ...,Hn) = det{M{L^)). 
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Example 3.2. Let K = Q{t), d = and consider the set of differential 
polynomials in K{xi, xz}{ui, U2}, 

U) = xi — 7t — ui — 3nii + 3u2 — tu2i 

F2{X, U) = X2 - Ui + Ui2 - 5U22 
F3(X, U) =Xs-Ui+ U12 - t^U2l. 

Then the set PS(Fi,F2,F3) contains L = 13 differential polynomials and 
9Res(-Fi, F2, -F3) = det(M(L)) where M{L) is the following L x L coefficient 
matrix of PS(Fi, F2, F3) . 



-t - 


-3 


-1 


-1 


























Xu 








-t 


-3 





-1 




















Xl3 














-t 


-3 
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-1 














X12 




















-t 


-3 
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-1 








Xu - 7 


























-t 


-3 


3 


-1 


xi — 7t 


-5 


1 











-1 




















X23 








-5 


1 











-1 














X22 














-5 


1 











-1 








X21 

















- 


-5 


1 











-1 


X2 





1 


-t^ 


- 


-6t 


-1 - 


-6 

















X33 











1 - 




- 


-At 


-1 


-2 











X32 

















1 - 


-t^ 





-2t 


-1 








X3I 























1 


-t^ 








-1 


X3 




be 


the 


square 


su 


bmatrix 


of 


size 


= 


10 


9/ M(L) 


obtained 



removing columns 1, 2, 13 and rows 1, 6, 10. Then dRes {Hi, H2, H3) = 
det (M(L^)). 

We show next some properties of 9Res(Fi, . . . , F„) and dRes'^{Hi, . . . , Hn) 
that will be used later in the paper. 

Proposition 3.3. // {Hi = 0, . . . , Hn = 0} has a nonzero solution then 
dRes'' {Hi,..., Hn) = 0. 

Proof. Let F be a differential field extension of K. Then every nonzero so- 
lution of {Hi = 0, . . . , Hn = 0} in F"~^ is a nonzero solution of the system 
{d^~°'~^Hi = 0, . . . , dHi = 0, Hi = \ i = 1, . . . , n}. If such a solution exists 
then the columns of M{L^) are linearly dependent on F. □ 

Remark 3.4. For n = 2, dRes^{Hi,H2) = and only if {Hi = 0,H2 = 
0} has a nonzero solution in ¥, a differential field extension 0/ IK (see f^, 
Theorem 3.1). Unfortunately, for n > 2 the condition dRes^{Hi, . . . , Hn) = 
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is not sufficient for the existence of nonzero solutions of the system {Hi = 
0, . . . , Hn = 0}. Let n = 3 and 

Hi{U) = Mil + M21, H2{U) =Ui+ U2, H'i{U) = Ml + Uu + M21. 

The first two columns of M{L^) are equal and therefore (9Res'^(ifi, iJ2) H^) = 0. 
The system {Hi = 0, H2 = 0, H^ = 0} has only the zero solution. 

We introduce some matrices that will be used in the proof of the next result 
and also in later results in the paper. Let S be the n x [n — 1) matrix whose 
entry {i,j) is the coefficient of Ujoi in Fi, i G {!,... ,n}, j G {1, . . . ,n — 1}. 
Let Si be the matrix obtained by removing the 2-th row of S. 

Remark 3.5. Note that the nonzero rows of the columns of M{L) ( resp. 
M{L^) ) corresponding to the coefficients ofujN (resp. ujn-i), j = ^, ■ ■ ■ — l 
are the rows of S . 

Proposition 3.6. Let M^.i be the L x (L — 1) principal submatrix of M{L). 
The following statements are equivalent: 

1. 9Res(Fi,...,F„) ^0. 

2. dRes^{Hi,...,Hn) ^ 0. 

3. rank(Mi_i) = L - 1. 

Proof. The matrix M{L^) is a square submatrix of M^^i so 3 implies 2. Given 
z G {1, . . . , n}, we call M\^_i the submatrix of M^^i obtained by removing the 
row corresponding to the coefficients of d^~°'Fi. Thus, if dR.esh{Hi, . . . , Hn) 7^ 
then by Remark 13.51 there exists k E {1, . . . ,n — 1} such that det(5'jt) 7^ 0. 
Furthermore, there exists a G N such that 

det(M^i) = {-lydReshiHi, . . . , det{Sk) + 0, 

which proves that 2 implies 3. 

Let K(X) be the quotient field of K{X}. The equivalence of 1 and 3 is 
obtained easily noting that the elements of the last column of M{L) are linearly 
independent on ]K(X) and so the last column of M{L) is linearly independent 
on K(X) from the first L — 1 columns. □ 

4 Implicitization of linear DPPEs by 
Carra'Ferros's differential resultant 

Let V{X, U), Fi, Hi be as in Section H Let D = K{X}, PS = PS(Fi, . . . , F„) 
and let ID be the implicit ideal of V{X,U). In this section, we use the 
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implicitization results in terms of characteristic sets given in [TT] to obtain 
implicit ization results in terms of differential resultants. We prove that if 
9Res(Fi, . . . , Fn) 7^ then 9Res(Fi, . . . , = is the implicit equation 



Recall that PS is a set of linear differential polynomials. Let [PS] be the 
differential ideal generated by PS, then it holds ID = [PS] n K{X} by pi] . 
Lemma 3.2. Let ^ be a characteristic set of [PS] and Aq = An K{X}. By 
Theorem 3.1, the implicit ideal is 



Furthermore, if | 1= 1 then Aq = {A{X)} where A is a characteristic 
polynomial of ID. By [6], Proposition 12, (9Res(Fi, . . . , Fn) G ID and it is our 
candidate to be a characteristic polynomial of ID. 

Let (PS) be the ideal in K[xi, . . . , Xijs-Oi.Uj, . . . , Uji^ \ i = 1, . . .n,j = 
1, . . .n — 1] generated by PS. To compute a characteristic set of ID we will 
use a Groebner basis Q of (PS) with respect to the ranking 71* on X UU that 
eliminates U with respect to X, that is d''x < d^*u, for all a; G X, u G f/ and 



Lemma 4.1. Let Q he the reduced Groebner basis o/ (PS) with respect to the 
ranking TZ* . 

1. Q = {Bq, Bi, . . . , Bl^i} where Bq < Bi < ■ ■ ■ < Bl_i with respect to the 
ranking TV and Qo = Q C\ K{X} is not empty. 

2. Let E{L) be the L x L matrix whose k-th row contains the coefficients 
of Bi^k, k = 1, . . . ,L as a polynomial in ©{f/}, and where the coeffi- 
cients are written in decreasing order with respect to the orderly ranking 
on U. Given the differential resultant matrix M{L) of Fi, . . . Fn, then 
det(M(L)) = {-lY det E{L) for some a G N. 

3. The cardinality of Qq is one if and only if 9Res'*(iJi, . . . , iJ„,) 7^ 0. 

Proof. Let M2L be the L x (2L) matrix whose fc-th row contains the coefficients 
of the /c-th polynomial in PS, as a polynomial in ]K{X U f/}, and where the 
coefficients are written in decreasing order w.r.t. TV . 



of P(X, U). 



ID = [PS] nK{X} = [^o]- 



k, k* G N| 



1 



qN-oi 



1 



M2L 



Ml 



-1 



1 



QN-o, 



a. 



•n 



1 



a. 



'n 
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Observe that we can find a minimal Groebner basis of (PS) performing 
gaussian elimination on the rows of M2L1 see [H], Section 7, exercise 10. In 
particular, the polynomials corresponding to the rows of the reduced echelon 
form E2L of M2L are the elements of the reduced Groebner basis Q of (PS). 

1. Observe that M2L has rank L. Therefore Q contains L elements < 
Bi < ■ ■ ■ < The sub matrix formed by the first L — 1 columns of 
M2L is Ml_i which has rank less than or equal to L — 1, therefore the 
set Qo = Q \~\ K{X} is not empty. 

2. The matrix E{L) is obtained by performing on M{L) the same row 
operations as the operations performed on M2L to obtain E2l- Therefore 
det(M(L)) = (-l)MetE(L) for some a e N. 

3. The cardinality of Qq is 1 if and only if rank(M£,_i) = L — 1. Equivalently, 
dRes^{Hi, . . . , i7„) 7^ by Proposition [SH 

□ 

To compute a characteristic set of [PS] we apply the algorithm given in [1], 
Theorem 6, that we briefly include below for completion. Given P G ]K{XUU}, 
the lead of P is the highest derivative present in P w.r.t. TZ*, we denote it 
by lead(P). Given P,Q ^ K{X U U} we denote by prem(P, Q) the pseudo- 
remainder of P with respect to Q, [15], page 7. Given a chain A = {Ai, . . . , At} 
of elements oiK{XUU} then prem(P, A) = prem(prem(P, At), {Ai, . . . , A-i}) 
and prem(P, 0) = P. 

Algorithm 4.2. Given the set of polynomials PS the next algorithm returns a 
characteristic set of [PS] . 

1. Compute the reduced Groebner basis Q of (PS) with respect to TZ* . 

2. Assume that the elements of Q are arranged in increasing order Bq < 
Bi < ■ ■■ < Bl-1 w.r.t. TZ*. Let A = {-Bo}- For i from 1 to L - 1 do, if 
lead(Pi) 7^1ead(Pi_i) then A := ^ U {prem(Pi, ^)}. 

Theorem 4.3. Given a system V{X, U) of linear DPPEs with implicit ideal 
ID. If dIles^{Hi, . . . , Hn) 7^ then ID has dimension n — 1 and 

(9Res(Pi,...,P„)(X) = 

is the implicit equation ofV{X, U). 



12 



Proof. Let Q be the reduced Groebner basis of (PS) with respect to the ranking 
TZ* and let i?o < -Bi < ■ ■ ■ < -Bl-i be the elements of Q. By Lemma 14.11 we 
have Qo = Q r\ K{X} = {-Bo} and so there exists a characteristic set A of [PS] 
such that ^0 = {-Bo}- Consequently, the dimension of ID is n — 1. 
By Lemma 14. II and the definition of the differential resultant, 

9Res(Fi, ...,Fn)= det(M(L)) = (-1)" det(E(L)), 

for some a G N, therefore (9Res(-Fi, . . . , F„) = cBq with c G K. This proves 
that 9Res(-Fi, . . . , F„) is a characteristic polynomial of ID and therefore the 
implicit equation of V{X, U) is aRes(Fi, . . . , F„) (X) = 0. □ 

5 Linear Complete Differential Resultants 

Let D be a differential integral domain, and let /j G D{f/} be a linear ordinary 
differential polynomial of order Oj, i = 1, . . . , n. For each jG{l,...n — l}we 
define the positive integer 

7j(/i, ■■■Jn)= min{oi - ord(/„ Uj) \ i e {1, . . . , n}}. 

Observe that, < 7j(/i, . . . , /„) < Oj for each z G {1, . . . , n}. We also define 

n-l 

7(/l, • • • , /n.) = X] 7i(/l, • • • , fn) 

i=i 

and it is easily proved that < 7(/i, . . . , fn) N — Oi, for alH G {1, . . . , n}. 

Observe that, when 7(/i, . . . , fn) 7^ then the set PS(/i, . . . , /„) defined in 
Section [31 is a set of L differential polynomials in L — 7(/i, . . . , /„) — 1 differ- 
ential variables. Then c}Res(/i, . . . , fn) is the Macaulay's resultant of a set of 
polynomials which are not complete in all its variables, thus c?Res(/i, . . . , /„) = 
0. If 7(/i,...,/„) = we will say that the set of differential polynomials 
{/i, . . . , /„} is complete. We will call 7(/i, . . . , fn) the completeness index of 

{/l; • • • ; fn}- 

We define next a differential resultant that generalizes Carra'Ferro's dif- 
ferential resultant in the linear case. Observe that the generalized differential 
resultant in deals with a different aspect. 

Definition 5.1. The linear complete differential resultant 9CRes(/i, . . . , /„), 
of n linear ordinary differential polynomials /i, . . . , /„ in n — 1 differential 
variables Mi, . . . ,Mn-i, is defined as the Macaulay's algebraic resultant of the 
differential polynomial set 

PS7(/i, . . . , fn) '■ = 

n 

{a^-°-Vi,...,5/.,/iM = l,...,^, N = J2oi, 7 = 7(/i, •••,/«)}. 

1=1 



13 



The set PS^(/i, . . . , /„) contains = ^"^^(A^ — Oj — 7 + 1) polynomials 
in the following set V of — 1 differential variables 

V = {uj, UjN-^j-j I 7i = Ijifi, ■■■,fn),j = ^,---,n-l}. (3) 

Observe that EJ=i(^ - 7^ - 7 + 1) = (n - 1)A^ - ^7 + n - 1 = - 1. 

Let hi G D{f/} be a linear ordinary differential homogeneous polynomial 
of order Oi, i = 1, . . . ,n. 

Definition 5.2. We define the linear complete differential homogenous resultant 

9CRes'*(/;.i, . . . , hn), of n linear ordinary differential homogeneous polynomi- 
als hi, ... ^h^ in n — 1 variables, as the Macaulay's algebraic resultant of the 
differential polynomial set 

PS;(/ii,...,M := 

n 

{d''-"^-^-'h„...,dhi,h,\z = l,...,n, N = Y,o^, 7 = 7(^i,---,M}- 

i=l 

The set PS!^(/ii, . . . , hn) contains = ^^^i{N — Oj — 7) polynomials in 
the following set of — 1 differential variables 

= {Uj, . . .,UjN-yj-y-l I 7i = lj{hl, ■ ■ ■,hn),j = 1 , . . . , W - 1 } . 

Observe that for 7(/i, . . . , fn) = (resp. j{hi, . . . , hn) = 0) it holds 

(9CRes(/i, ...,/„) = (9Res(/i, ...,/„), 
denes'" {hi, ...,hn) = dRes''{hi, . . . , /i„). 

We introduce next the matrices that will allow the use of determinants 
to compute 9CRes(/i, . . . , /„) and c?CRes'^(/ii, . . . , /i„). For i = 1, . . . , n, 7 = 
lifi, ■■■,fn) (resp. 7 = 7(/ii, . . . , hn)) and /c = 1, . . . , iV - Oj - 7 (resp. k = 
1, . . . , N — Oi — •y — 1) define de positive integers, 

l{i,k) = (i-l)(iV-7)-^0i + i + /c, 

h=l 

(f'li^k) = {i - 1){N - - 1) -J^Oi + i + k. 

h=l 

Then l{i,k) G {!,..., L^} and l'^{i,k) G {!,..., L^J}, with the appropriate 
value of 7 in each case. 

Let M{L^) be the x (resp. x L^) matrix containing the coefficients 
of a^-°-T-^7i ( resp. d^~°^~^-''-^hi) as a polynomial in D[V] (resp. in ©[V^]) 
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in the k)-th row (resp. k)-th row), where the coefficients are written in 
decreasing order with respect to the orderly ranking on U. Hence, M{L^) (resp. 
M{L^)) is a matrix over D that we call the complete differential (homogeneous) 
resultant matrix of /i, . . . , /„ (resp. hi, ... , hn). In this situation: 

9CRes(/i,...,/„) =det(M(L^)), 
dCRes\hi, ...,K)= det{M{L^)). 

The next matrices will be used in the proofs of the results in this section 
and also in later results in the paper. Let be the n x (n — 1) matrix whose 
entry is the coefficient of Mjoi-7j /«; ^ {I5 • • • > ''T'}: 3 ^ {1; • • • > ''^ ~ 1} 
For i G {1, . . . n} let S.^i be the (?i — 1) x (ri — 1) matrix obtained removing the 
i-th row of S.y. 

Remark 5.3. Let us suppose that hi is the homogeneous part of fi, i = 
1, . . . , n, then 7 = 7(/i, . . . , /„) = 7(^1, • • • , hn)- Note that the nonzero rows 
of the columns of M{L^) ( resp. M{V^)) corresponding to the coefficients of 
UjN-'^ (resp. UjN-'y-i), j = ^, ■ ■ ■ ,n — 1 are the rows of Sj. 

Theorem 5.4. 1. If the system {hi = 0, . . . , hn = 0} has a nonzero solu- 
tion then (9CRes'^(/ii, . . . , hn) = 0. 

2. Let us suppose that hi is the homogeneous part of fi, i = 1, . . . ,n. Let 
Ml^_i be the x (L^ — 1) principal submatrix of M{L^). Then, 

(a) dCRes^{hi, . . . , hn) 7^ if and only z/rank(Mi^_i) = — 1, 

(b) if 9CRes(/i, ...,/„) ^ then dCRes^{hi, . . . , /i„) ^ 0. 

Proof. 1. The proof is analogous to the proof of Proposition 13. 3[ 

2. The matrix M{L^) is a square submatrix of Mi_^_^ . Given i G {1, . . . , n}, 
we call M|^_^ the submatrix of Mx,^_i obtained by removing the row cor- 
responding to the coefficients of d^~°^~"' fi. Thus, if dCKes^{hi, . . . , hn) 7^ 
then by Remark I3.5l there exists /cG{l,...,n— 1} such that det(S'^fc) 7^ 
0. Furthermore, there exists a G N such that 

det(M^^„,) = (-l)«9CRes'^(/ii, . . . , /i„) det{S^k) 7^ 0, 

which proves the equivalence (a). Since 9CRes(/i, . . . , /„) 7^ imphes 
rank(Mi;,^_i) = — 1, statement (b) follows. 

□ 

Corollary 5.5. Let Fi and Hi be as in Section [1 then (9CRes(Fi, . . . , Fn) = 
if and only if dCRes''{Hi, ...,Hn)=0. 
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Proof. By Theorem \EM 2(a), 9CRes(Fi, . . . , F„) ^ if and only if it holds 
rank(Mi^_i) = — 1. The proof is analogous to the proof of Proposition 
1331 □ 

Example 5.6. Let K = Q(t), d = and consider the set of differential 
polynomials in K{xi, X2, x^}{ui, ^2}, 

U) = Xi - t - Ui - U2 - 2U22 

F2(X, U) = X2-t^ - 2ui - U2- U22 
F^i^X, U) = X3 - 5 - ui - 3uu - U21 - U23. 

Then the set PS(Fi,F2,F3) contains L = 17 differential polynomials and 
9Res(-Fi, F2, -F3) = det(M(L)) = because the columns of M{L) corresponding 
to the coefficients of Un and Uiq are columns of zeros. 

We have 71 = 2 and 72 = so 7 = 2. Let us compute (9CRes(Fi, F2, F3). 
The set PS^(Fi, F2, F3) contains L^ = 11 differential polynomials in the dif- 
ferential variables V = {^25, U24, U23, M13, U22, M12, M21, Uu, U2, Ui} written in de- 
creasing order w.r.t. the ranking on U . 



M{L^) = 



-2 





-1 


-1 

























-2 








-1 


-1 














X12 








-2 











-1 


-1 








Xu - 1 














-2 











-1 


-1 


Xi — t 


-1 





-1 


-2 




















X23 





-1 








-1 


-2 














X22 - 2 








-1 











-1 


-2 








X21 - 2 1 














-1 











-1 


-2 


J.2 
X2 - t 


-1 





-1 


-3 





-1 














X32 





-1 








-1 


-3 





-1 








X3I 








-1 











-1 


-3 





-1 


X3-5 



The matrix M{L^) is obtained by removing rows 1, 5, 9 and columns 1, 4, 
11 of M{L^). Then, 

aCRes(Fi, F2, F3) =4x3 - 8 - 8x21 + I2t + 12x32 + 4xi3 + 4xii - 20x23 
- 8x22 + 4xi2 + 4xi - 4x2 + 
dCRes^{Hi,H2,Hs) = - 4. 

6 Implicitization of linear DPPEs by linear 
Complete Differential Resultants 

Let V{X,U), Fi, Hi be as in Section [2] and let ID be the implicit ideal of 
V{X,U). In this section, we prove the main result of this paper, namely if 
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9CRes {Hi, . . . , Hn) 7^ then 5CRes(Fi, . . . , F„)(X) is the imphcit equation 
of V{X, U) and hence the results in Section H] are extended. 

Let 7 = 7(Fi, . . . , F„) = -f{Hi, ...,Hn) and V as in ([3]). Let PS^ = 
PS^(Fi, . . . D = K{X} and let (PS^) be the ideal generated by PS^ in 
K[xi, . . . , XiN-Oi-'y, Uj, ■ ■ ■ , UjN-'y \ i = I, . . . H, j = 1, . . . Ti — 1]. Let 71* be the 
ranking on X U U defined in Section |H 

Lemma 6.1. Let Q he the reduced Groebner basis of (PS-y) with respect to the 
ranking 71* . 

1. Q = {Bq, Bi, . . . , Bl^^i} where Bq < Bi < ■ ■ ■ < Bl^^i with respect to 
the ranking 7Z* and Qo = Q r\ K{X} is not empty. 

2. Let E{L^) be the x matrix whose k-th row contains the coefficients 
of Bi^-k, k = 1, . . . , as a polynomial in D[V], and where the coeffi- 
cients are written in decreasing order with respect to the orderly ranking 
on U . Given the differential resultant matrix M[L^) of Fi, . . . Fn, then 
det(M(L^)) = (-1)" det E{L^) for some aeN. 

3. The cardinality of Qq is one if and only if dCRes''{Hi, ...,Hn)^0. 

4- If dCRes'^{Hi, . . . , Hn) 7^ then A = {Bq, Bi, . . . , B^-i} is a character- 
istic set of [PS^] and Ao = AU K{X} = {Bq}. 

Proof. The proof of 1,2 and 3 is analogous to the proof of Lemma 14. 1[ 

If dCRes^{Hi, ...,Hn)^0 then E{L^) is an upper triangular matrix of 
rank whose — 1 principal submatrix is the identity. Following Algorithm 
14.21 we obtain the characteristic set A = {Bq, Bi, . . . , Bn-i} of [PS^]. Since 
go = Gn K{X} = {Bo} then Ao = {Bo}. □ 

The next theorem gives an explicit formula of (9CRes(Fi, . . . , Fn) in terms 
of dCKes^{Hi, . . . , Hn). Given i G {1, . . . , n}, let S^i be the matrix defined in 
Section [5l 

Theorem 6.2. 1. There exists P G ID and k E {1, . . . ,n} such that 

( a) ord(P, Xk) = N - Ok - 7 

(b) ord(P, Xi) < N ~ Oi — -f , i = 1, . . . ,n. 

2. For any differential polynomial P as in statement 1 it holds 

aCRes(Pi, ...,Fn) = - det(5^fc)aCRes^'(ifi, . . . , Hn)P{X). 

a 

with a = (-1)'^^^^^ , aeN. 



17 



Proof. Given the reduced Groebner basis Q of (PS^) w.r.t. TZ* let us take a 
polynomial B & Qq = Q n K{X}. Let /? = min{A^ — — 7 — ord(i?, Xi) \ i = 
1, . . . , n} = iV — Ofc — 7 — ord(i?, Xk) for some /c G {1, . . . , n}. Then P = d^B 
verifies the conditions in statement 1. 

Since P e ID = [PS^] n K{X} there exist J^i, . . . , J^„, G such that 
P(X) = J'^{F^{X, U)) + ■ ■ ■+J'n{Fn{X, U)). Then P{X) = ^i(Ti(X)) + • • ■ + 
TniTn{X)) and Ti{Hi{U)) + ■ ■ • + JF„(if„(f/)) = 0. As a consequence, we can 
perform row operations on M{L^) to obtain a matrix of the kind 



0...0 

0...0 
0...0 



0...0 P{X) 

* ... * * 



* ... * 



To be more precise, we reorder the rows of M{L^) so that the coefficients of 
g^-Ok-if^ are in the first row, rows 2 to n are the rows containing the entries 
of the matrix S^k and rows + 1 to are the rows containing the entries 
of M{L^). Then multiply the first row of the obtained matrix by -q^ 



dP 



'^k N — oj^ —7 



Finally, replace the first row by the coefficients of P{X) as a polynomial in 
D{t/} written in decreasing order w.r.t. the ranking on U, that is, all zeros 
and the last entry equal to P{X). 

Let a = (— , for some a G N then 



dP 

k N-Of.-j '■ 



adet(M(L^)) = det(5fc)(9Res!;(ifi, . . . ,if„)P(X). 



□ 



Theorem 6.3. Given a system V{X, U) of linear DPPEs with implicit ideal 



ID. If dCRes" {Hi, 



Hn) 7^ then ID has dimension n — 1 and 
aCRes(Pi,...,P„)(X) = 



is its implicit equation. 

Proof. By Lemma 16.1( 4). a characteristic set of ID is ^0 = {-^o}- Then the 
dimension of ID is n — 1. By Lemma 16.1( 2). and the definition of the linear 
complete differential resultant. 



aCRes(Pi 



det(M(L^)) = (-l)"det(E(L^)), 



for some a G N. Therefore 9CRes(Pi, . . . , P„) = cBq with c G IK. This proves 
that c}CRes(Pi, . . . , P„) is a characteristic polynomial of ID and therefore the 
imphcit equation of U) is aCRes(Pi, . . . , Fn){X) = 0. □ 
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7 Some results on properness 

Let V{X,U), Hi, Pi and Cij be as in Section [2l In this section, we give 
some results related with the inversion problem. To be more precise, we study 
conditions on the differential operators Cij so that V{X, U) is a set of proper 
DPPEs. 

We gather below some definitions that will be needed in this section and 
that were used in [TT], Section 6, to study the inversion problem in terms of 
characteristic sets. The image of 'P{X, U) is the set 

IM = {(r/i, . . . , r/„) e I 3(ri, . . . , r„„i) G E^-^ with r,, = P,{n, r„_i)}. 

The inversion problem says: given (771, . . . , ?7„) G IM, find (ri, . . . , r„_i) G E"^^ 
such that rji = Pi(ri, . . . , r„_i). We call inversion maps for V{X, U) to a set of 
funtions g^, . . . , g^-i in {X} such that 

Uj = gj{xi, . . . j = 1, . . . ,n - 1. 

A set V{X, U) of DPPEs is proper if for a generic zero (Fi, . . . , F^) (and hence 
most of the points) of the implicit variety Zero(ID), there exists only one 
(ti, . . . , r„_i) G E"~^ such that Fj = Pj(ri, . . . , r„_i). By [TT], Lemma 3.1, 
(Pi(?7), . . . , PniU)) is a generic zero of the implicit ideal ID. 

Proposition 7.1. Let as suppose that (9CRes^ (i/i, . . . , if„) 7^ 0, then the next 
statements hold. 

1. The set V{X, U) of DPPEs is proper. 

2. Furthermore, there exist a set of linear inversion maps for V{X, U), 

f/i,...,f/„„iGK{X}. 



Proof. By Lemma 16.11 then E{L^) is an upper triangular matrix of rank 
whose L^ — 1 principal submatrix is the identity. Then, A = {Bq, Bi, . . . , 
is a characteristic set of [PS^] where 

B,{X, U) = u, - U,{X), J = 1, . . . , n - 1, 

for linear differential polynomials Uj G K{X}. By [11], Theorem 6.1, the set 
V{X, U) of DPPEs is proper and Ui, . . . , Un-i is a set of inversion maps of 

v\x,u). □ 

If IK is not a field of constants with respect to 5, then is not commu- 
tative but 

dk — kd = d{k) 
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for all A; G K. The ring ]K[d] of differential operators with coefficients in K is 
right euclidean (and also left euclidean). Given C,C G K[d], by applying the 
right division algorithm we obtain g, r G IK[9], the right quotient and the right 
reminder of C and C respectively, such that C = qC'+r where deg(r) < deg(g). 

If Cj G K[d] is the greatest common right divisor of Cij, . . . , Cnj then there 
exists C'ij G K[d] such that Cij = C^jCj, i = 1, . . . ,n. Then we write 

Cj = gcrd(£ij, . . . , Cnj), j = 1, . . . , n - 1. 

By Remark 12. 2^ then Cj ^ 0. If Cj G K, then we say that C\j, . . . , C^j are 
coprime and we write 

{^C\j, . . . , Cnj) = 1- 

Theorem 7.2. A necessary condition for the set V{X, U) of DPPEs to he 
proper is 

{Cij, . . .,Cnj) = 1, j = 1, . . . ,n - 1. 

Proof. Let us suppose that there exists k G {1, . . . ,n — 1} such that Ck = 
gcrd(£ifc, . . . , Cnk) is non constant. Then there exists a nonzero element r/ G K 
such that Cikirj) = 0, z = l,...,n. Define the element 



n-l 



U + r] = ...,Uk + ri,..., Un-i) G E' 

Recall that Pi{U) = ai - Y^'^Zl ^ij{uj), then {Pi{U + r/), . . . , P„(t/ + rj)) = 
(PiiU), Pn{U)). Thus by definition, V{X, U) is not proper. □ 



Remark 7.3. For n = 2 the condition in Theorem \7.^ is also sufficient but 

df 



for n > 3 this is not true. Let K = C(t) and d = Let us consider de system 



of linear DPPEs 




where 



2Ui + Mil +U2 + U22 

Ul + nil + Ui2 + U2 — U22 

ui + 2nii + ^2 + U2i 



Cu = 2 + d C2i = l + d + d^ £31 = 1 + 29 
Ci2 = l + d £22 = 1 + 92 £32 = 1 + d. 



We compute a characteristic set A of [PS] with Algorithm \4.2\ 

A = {Xi2 - 2X22 - 2X21 -X2+ X33 + X32 + X31 + X3, 

U2 + 3/2ni + xii - l/2xi - 2x2i - X2 + X32 + l/2x3i + 1/2x3, 

Uu - Ul- Xu +Xi + 2X21 - X32 - Xs} 

Then by /7i]/ . Theorem 6.1, the system is not proper but {Cij , C2j , C^j) = 1, 
J = 1,2. 
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We define a new system of DPPEs liaving tlie same implicit ideal than 



where H'^iU) = YT'l C-'.,{u,) and i^'(f/) = T,(X) + H[{U). 

Proposition 7.4. Let ID' he the implicit ideal of the setV'{X,U) of DPPEs. 
Then ID = ID'. 

Proof Given / E ID' then f{P[{U), . . .,Pn{U)) = 0. In particular, for r] = 
. . . , /:„_iK_i)) then = fiPiiv), P'M) = fiPiiU), Pn{U)). 
Therefore ID' C ID. 

Now let us suppose that / G ID = [PS] fl ]K{X}. Then, there exists jFj e 
i = such that /(X) = J-i(Fi(X, + - ■ ■ + [/)). Thus, 

+ ■ ■ ■ + J'n^HniU)) = and /(X) = J^iiT^^X)) + ■■■ + .F„(T„(X)). 
As a consequence, for each j G {1, . . . , n — 1} we have J^i{Cij{uj)) + ■ ■ ■ + 
J^n{Cnj{ui)) = • Thus {TiC[j + ■■■ + TnC^nj)l^o = and by Remark O 
Cj 7^ so the differential operator + ■ ■ ■ + TnC-'nj = 0. We conclude that 

Ti{H[{U)) + ■ • ■+J^n{K{U)) = which implies f{Pi{U), Pl^{U)) = and 
proves ID C ID'. □ 

The following corollary follows directly from Theorem I6.3l and Proposition 



Corollary 7.5. Given a system V{X, U) of linear DPPEs with implicit ideal 
ID. // 9CRes'^(if(, . . . , H'^) 7^ then ID has dimension n — 1 and 



is its implicit equation. 

Example 7.6. The method to compute the implicit equation of the system 
V{X, U) given by the equations xi = U1+U11 + M2 + M21; X2 = ^(mii + M12) + ^22; 
3^3 = U1 + U11 + U21 would be the following. First, compute the system V {X , U) 
as in (jl]), which is the system in Example Then the implicit equation of 
V{X, U) IS 9CRes(Fi', F2', F^)(X) = {t - l)xi2 - tegi - {t - l)x32 + X2 = 0. 

8 Treatment of special cases 

We give an explicit expression of the implicit equation of 'P(X, U) in terms 
of the differential operators defining the DPPEs for n = 2 and with some 
restrictions for n = 3. 



V{X,U) 




El 



aCRes(F;,...,F^)(X) = 
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8.1 Case n=2 

Given differential operators >Ci, £2 G ^[d] we consider the set of DPPEs 

Let Hi{u) = Ci{u), H2{u) = C2{u) and Fi(xi,X2,u) = xi — ai + Hi{u), 

F2(xi,X2, u) = X2 - a2 + H2{u). 

Remark 8.1. By JB^, Theorem 2, it holds dRes'^{Hi, H2) ^ if and only if 

If K is a field of constants with respect to 9, that is d{k) = for all A; G K 
( for example d = and IK = C) then C1C2 — C2C1 = 0, i.e. the differential 
operators commute. If IK is not a field of constants with respect to d, then 
K[9] is not commutative. We need the next commutativity result to prove the 
results in this section. Most likely it was studied by Tsirulik in [19] but we 
have not been able to find this paper. We give an algorithmic proof. 

Lemma 8.2. Let Ci, C2 G K.[d] such that {Ci, C2) = 1. Then there exist 
T^iiT^2 G with deg{T>i) < deg{Ci) — 1, i = 1,2 such that 

(£2 - I?2) A - (A - A = 0. 

Proof. If Ci and C2 commute then Vi = 0, i = 1, 2. Let us suppose that 
C1C2 - C2C1 ^ 0. Let £1 = E°io'/'i^' and £2 = EjLo^i^' ^i^h e K. 

By equation (1.2) in [2], we have 

01+02 

£i£2 - £2£i = ^ (cfc - Ck)d^ where 

fc=0 

min{oi,k) oi min(o2,k) 02 

Ck= Yl X]'^*^^'"'^^^'^-)' ^'^= E J2'^jd'-^~'\<l>k-s). 

s=max{0,k~O2) i=s s=max{0,k—oi) j=s 

Then Coj+oa = Coj+oa = 'Poi4'o2- Therefore the degree of the operator C1C2 — 
£2£i is less than or equal to Oi + 02 — 1. 

We want to find = E°io^ and V2 = J2T=o^ Pjd^ with at, pj E K 
such that ■Pi£2 — ^^2£i = £i£2 — £2£i- We have 

'D1C2 - 'D2C1 = ^ (7fc - 'lk)d'' where 

k=0 

min{oi~l,k) oi — 1 min{o2 — l,k) 02 — 1 

s=max{0,k—O2) s=max(0,fc— 01) i=s 
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Let us consider the system of Oi + 02 equations 

7k - % = Ck - Ck, = 0, . . . , oi + 02 - 1 (6) 

in the N = oi + 02 unknowns a^, i = 0, . . . , Oi and Pj, j = 0, . . . , 02. The 
coefficient matrix of the system ([6]) is the differential homogeneous resultant 
matrix M{L^) of Hi, H2, with = N. Thus the system has a unique solution 
since det(M(L'')) = dRes^{Hi, H2) ^ 0, by Remmark [HIl □ 

Theorem 8.3. Given a system V2{xi,X2,u) with implicit ideal ID. 

1. The following statements are equivalent. 

(a) 'P2{xi,X2,u) is proper. 

(h) {Ci,C2)=l. 

(c) dRes''{Hi,H2) ^ 0. 

2. If{Ci, C2) = 1 then the dimension of ID is 1 and dIles{Fi, F2){xi, X2) = 
is the implicit equation of V2{xi, X2, u). Furthermore, there exist Vi, V2 G 
K[9] such that 

dRes{Fi,F2){xi,X2) = 

{-irdRes\Hi, H2)[{C2 - V2){xi - a^) - {Ci - I?i)(x2 - 02)], 
for some a G N. 

Proof. 1. It follows from Proposition 17.11 Theorem 17.21 and Remark 18.11 

2. Observe that for n = 2, 7(^1,^2) = j{Hi,H2) = 0. By statement 
1 and Theorem 16.31 dRes{Fi, F2){xi,X2) = is the implicit equation 
of V2{xi,X2,u). By Lemma 18.21 there exist Vi,V2 G ]K[d] such that 
(£2 - T>2)Ciiu) - (£1 - Vi)C2iu) = 0. Let 

P(xi, X2) = {C2 - V2){xi - ai) - (£1 - Vi){x2 - 02). 

Since dRes^{Hi, H2) 7^ then there exists k G {1, 2} such that det(S'fc) = 
dP/dxkN-o,, 7^ 0. The result follows from Theorem 16. 2[ 

□ 

Algorithm 8.4. Given the system V2{xi,X2,u) of DPPEs this algorithm re- 
turns its implicit equation A{X). 

1. C := gcrd(£i,£2)- 

2. If C then compute £■ such that Ci = C'^C and set Li := z = 1, 2. 
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3. Compute T>i and T>2- 



4. A{X) := (£2 - V2){xi - ai) - (£1 - I?i)(x2 - aa). 

The Maple packages OreTools, DEtools and Ore-algebra can be used to 
compute gcrd(£i,£2), [I]- Also, we can compute Vi and V2 following the 
algorithm in the proof of Lemma 18.21 

8.2 Case n=3 

We restrict to the case K = C and consider the system of DPPEs 

{xi = ai- £ii(ni) - £i2(m2) 
X2 = a2-/:2i(ni)- £22(^2) (7) 
X3 = as - C3i{ui) - £32(^^2) 

where X = {xi,X2,X3}, U = {^1,^2} and Cij G C[d]. These differential 
operators commute, therefore the polynomial 

P{X) =£2i£32(xi - ai) - £22£3i(a;i - ai) - £ii£32(x2 - 02)+ 

+ £i2£3l(x2 - 02) + £ll£22(a;3 - 0.3) - Ci2C2l{x3 - 03) 

belongs to the implicit ideal of Vsi^X, U). 

For i = 1, 2, 3, we have U) = ai+Hi{U) and Hi{U) = £a(wi)+£i2(^^2) 

of order Oj. Let 7 = j{Fi, F2, F3) and let S^i be the matrix defined in Section 
[51 It is easily obtained that 

dP 

det(^70 = <T- • 

If det(S'^j) ^ for some i G {1, 2, 3} then it follows from Theorem 16.21 that 
aCRes(Fi,F2,F3)(X) = {-irdCRes\H^, H2, H3)P{X), 

a G N. 
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